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ABSTRACT 



We study generalized matrix models corresponding to n-point Virasoro conformal 
blocks on Riemann surfaces with arbitrary genus g. Upon AGT correspondence, these 
describe four dimensional J\f = 2 SU(2) n+39 ~ 3 gauge theories with generalized quiver 
diagrams. We obtain the generalized matrix models from the perturbative evaluation of 
the Liouville correlation functions and verify the consistency of the description with 
respect to degenerations of the Riemann surface. Moreover, we derive the Seiberg- 
Witten curve for the N = 2 gauge theory as the spectral curve of the generalized 
matrix model, thus providing a check of AGT correspondence at all genera. 
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1 Introduction 



The distinctive feature of M-theory is the description in geometrical terms of non-perturbative 
phases of superstrings. This approach is very effective for local geometries, where the dynamics 
of gravitational degrees of freedom gets decoupled and we gain a framework for the description 
of non-perturbative gauge theory dynamics. M-theory beautifully encodes the Seiberg-Witten ge- 
ometry of four dimensional N = 2 theories in terms of M5-brane compactifications El El El . In 
particular in [0 a full class of generalized quiver gauge theories has been described in terms of 
multiple M5-brane systems covering a generic punctured Riemann surface C 9jn . For example, for 
Co, n and Ci „ one recovers Witten's constructions of linear and circular quivers in the appropriate 
degeneration limits. 

In this context a very intriguing relation between the partition function of four dimensional 
SU(2) n+3a ~ 3 superconformal J\f — 2 gauge theories |4) and Liouville theory on C 9 n has been 
discovered in 0. This proposal has been a subject of intensive investigations and refinements 
from different viewpoints. Evidence for this conjecture as well as complete proofs for some cases 
can be found in j6l |7). Extensions to higher rank gauge groups and Toda field theories were 
introduced and discussed in |[8). The refinement of the correspondence in presence of gauge theory 
observables has been presented and studied in [|9l [10l QT). Moreover, some arguments for the 
derivation of the AGT correspondence were proposed in the M-theory context in lfT2l and via 
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matrix models in (13 [II [13 [IB El EC [IS IIS Ell E2|. 

Here we would like to address this correspondence from a complementary point of view, ex- 
plaining how to recover the geometry of the M-theory set-up and the Seiberg-Witten data in the 
wildest generality. To this end we derive a generalized matrix model from Liouville theory on C g ^ n 
and study its large N limit recovering the gauge theory Seiberg-Witten curve as its spectral curve. 
This provides a check of the AGT conjecture at all genera. 

In Section 2 we derive the generalized matrix model - as extended Selberg integrals - starting 
from the Coulomb gas representation of the residues of the perturbative Liouville theory correla- 
tors. The matrix model potential that we get has the form anticipated by Ifl3ll and in the elliptic 
case it coincides with the one derived in |[20l . 

In Section 3 we discuss the stability of this picture and its consistency with respect to the 
degeneration of the curve C 9i „ in general and present the degenerations of punctured tori as an 
explicative example. 

In Section 4 we analyze the large N limit and show how, by consistently adapting to our case 
the standard matrix model techniques, one gets a spectral curve in terms of quadratic differentials 
on C g .„ precisely reproducing the Seiberg-Witten curve and differential proposed in Q. 

We leave our concluding remarks to Section 5 and devote an Appendix to the detailed study of 
the degenerations of the C 2 ,o Seiberg-Witten data. 

2 From Liouville theory to generalized matrix model 

In this section we derive the generalized matrix model which corresponds to the n point conformal 
block on a Riemann surface C g of genus g. We derive it from the perturbative calculation of the 
correlation function of the Liouville theory by following the discussion in ||2~3~1 . 

The n-point function of the Liouville theory on C g is given by the following path integral 




Liouville on C g 




(2.1) 



where the Liouville action is given by 




d a <pd b <p + QR<p + 4TTfie 2b *). 



(2.2) 



We divide the Liouville field into the zero mode and the fluctuation (p(z, z) = O + <j>(z, z \ obtain 
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ing 

k=i 

§ = h I d2z Vg(9 ab 9a^ + QR4>). (2.4) 



where 



We can integrate out the zero mode O as 



+ "# e V " 0/<i2z#2i V 2 ( 9 - 1 )*e-^^ m ' = N) ( [ d 2 z^/ge 2 A N , (2.5) 



where we have used J d 2 Zy/gR = Anx = 87r(l _ and iV is defined as 



Therefore, the n-point function can be written as 

N-pr at\ r / r „\ N n 



Ie^+i^-^- (2 - 6) 



A = ^ r( N) I V^e-~ s ( I d 2 z^ge 2b ~A f[e~ 2m ^ w ^. (2.7) 



k=l 



When N E Z> , the correlator diverges due to the factor T(—N). The residues A N at these simple 
poles are computed then in perturbation theory in b around the free scalar field action (12.41) . From 
now on, our convention is that 



(•••>f 

which leads to 



Jv^^^ 2 ^" 9 ^..., (2.8) 



A N = L^- / e -£ J<*V^« f[ J d 2 Zl ^e 2b ^f[e- 2m ^ w A . (2.9) 

\ i=1 k=l I free on C g 

The condition (12.61) ensures momentum conservation in the free theory. 

Here we choose as a reference volume form d 2 z^/g = \co(z)dz\ 2 where u(z) is the coefficient 
of a reference holomorphic differential. This differential has 2g — 2 zeros, which we denote by 
(/ = 1, • • • 2g — 2). Then, the first factor in the expectation value of (12.91 ) becomes 

or or 2g ~ 2 2g ~ 2 

f - / d 2 z^z)dd\og\uj\ 2 ^f- / d 2 z^z)^W5 2 (z-^i) = Q ( 2 - 10 > 
71 J 71 J 1=1 1=1 
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where we have used R = —(2/ v fg)dd\og y /g. Thus, we obtain 



.4 



A 1 



-/*) 



AT /2g-2 



2bN\ 



AT 



e 

J=l " i=1 fc=l / free on C 



J e -W(«*) 



(2.11) 



The i'-point function of the free theory on C g is given in the factorized form as Il24ll25ll26 



(detlmr) 1/2 6(J^h 



i=i 

oo 9 



free on C 



oo 



a=l 



Zi.zy)* 4 *^ 2 exp 



1=1 



27T? ^PaPbTab + 2^^^^ 

a, 6 a,? 



(2.12) 



where r ab is the period matrix, E [z^ Zj) is the prime form, {co a } is a basis of normalized holomor- 
phic one-forms, and p a is interpreted as the momentum flowing through the a-th A-cycle. 

Using the explicit expression (12.121) for (12.91) . we find that the residue An of the n-point func- 
tion of the Liouville theory reduces to the following integral 

C (w, m, £, p) exp I 2ni ^ p a p b T ab 

V 0,6 

exp I 4Trb^2p a / u a 



9 p+00 
A N oc [ / dp a 

N 

n / d*zMzi)\ 2+2b2 



i=l 



\[E{z u z 3 )- 2b2 \[E[ 



:2bmi 



i, 2 



(2.13) 



i<j i,k 

where the factor C(w, m, £, p) was defined as 

n 2 

C(w,m,£,p) = [[uj(w k )^T [[E(w kl w £ ) £ [[E(^,w k ) 

I k=l k<l I,k 



(2.14) 



\E(0,Or^ exp 



As in the torus case 112011 . it is not straightforward to factorize the integrals over the Riemann 
surface into holomorphic and anti-holomorphic integrals for generic N. However this is easily 
performed in the large N limit. Indeed, the last two-lines of (12.131) can be written as 



ni2 I — W\2 
a Zi l/zess I ~ 



\dz. 



-w 



(2.15) 
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where /i and W are 



W 



■1-6* 



A 1 



i=i \fc=i 



J] ^(^i, Zj)~ 2b2 Y[ E(zi, z *) 2b ^ kmk ^ 9s , (2.16) 

l<i<j<N i 

W«) , (2-17) 

a=l J ) 



where we have chosen a base point z* in order to split the measure from the potential and we have 
rescaled the parameters as m k — > m k /g s and p a — > p a /g s - Notice that the term in the square 
brackets in (12.161) is independent on the zeroes of the conformal factor, ensuring therefore that the 
generalized matrix model correctly embodies the conformal symmetry of Liouville theory. 

The large N limit amounts to take g s — > keeping g s N, b, m k and p a finite. In this limit, the 
conditions for the criticality are given by 



E'( 



k=i 



E\Zi . Zj t 



dZi - ) m k — -dZi - 2n > p a u a {Zi) = 



(2.18) 



a=l 



where E'(zi,z 2 ) = d zl E{z\, z 2 ). The conditions obtained from the ^-derivatives are just the 
complex conjugate of (12.181) . It is remarkable that the conditions for criticality are separated into 
holomorphic and anti-holomorphic equations, which implies that the integrals over the Riemann 
surface in (|2.13l) can be factorized into holomorphic and anti-holomorphic integrals in the large N 
limit as stated in (12.151) . We are therefore left with the following matrix-like integral 



N 



w,m,p,v) = 



\dzi 



8=1 



iV 



i,i 



-i-b 2 



I b s-^ / v^ , E{Zi,w k ) 
ys i=\ \k=i 



E( Zi ,z* 



n 

l<i<j<N 
9 

f47T^J> 



E(Zi, Zj) 



-2b 2 



ri £ < 



,*\ 2b 2~2 k m k /g a 



0=1 



(2.19) 



where w = {w k } , m = {m fc }, p = {p a } and v = {u a } are the filling fractions v a = bg s N a 
which specify the holomorphic integral above. The integrand in (12.191) is a proper one-from in 
each variable z$ on the covering space of the Riemann surface due to momentum conservation 
(12.61) . The matrix model potential that we find is in the form anticipated by [fT3l . 



In order to count the number of moduli of our matrix model we should note that there are 
n + 2g — 3 independent filling fractions: naively the number of critical points of the action is 
2^ — 2 + 71+1. However, there are constraints coming from the fact that we are free to move 
the base point z* and that we have specified the residue at the base point as above by using the 
momentum conservation. The latter is equivalent to the constraint on the sum of filling fractions 
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v a = bg s N. These constraints reduce the number of moduli by two, thus giving the correct 
counting. The paths of the integrals are defined such that only the solution of (|2.18l) labeled by 
the fixed filling fractions {u a } contributes to the integrals. The measure factor in (12.191 ) can be 
regarded as a generalization of the Vandermonde determinant. The differential dzd z W has simple 
poles with residues ({2m k }, — 2 J2k m k) at the points ({wk},z*). 

The integral in (12.131) is then obtained by integrating (12.191) and its complex conjugate over the 
filling fractions. Thus, in the large N limit, An becomes 



n+2g-3 



/oo a p 1 "a " 
\dp a I dv k 
-oo „_i J i i 



a=l " fe=l 



9 9 

exp ( 2iri^2^2p a r ab p b ] <7(w,m,p,£)Z^ s,n (w,m,p, v ) 

a=l 6=1 



(2.20) 



At the level of the generalized matrix model, the filling fractions v a are free parameters. Together 
with p a (a = 1 , . . . , g) which are independent parameters in the potential, we have totally n + 3g — 3 
independent moduli which are identified with the internal momenta in the Liouville conformal 
block and then with the Coulomb moduli of the gauge theory. The explicit correspondence of these 
parameters with the internal momenta is discussed in detail in section l4~Tl Under this identification, 
we see from C£3& that Zjf' n (v, m, p, v) is proportional to the conformal block of Liouville theory 

ma. 



3 Degenerations 

To study the behavior of the generalized matrix model when approaching perturbative corners 
in the space of gauge couplings, we have to study what happens when we degenerate C 9jn . The 
degeneration is usually described by using the plumbing fixture decomposition of the curve. Let 
lA t be the annulus 

Ut = {(z,w)\zw = t; \t\ < \z\ < 1; |t| < \w\ < 1} 

which as t ~ describes the squeezed cylinder. The curve undergoes the decomposition C g n = 
C g -i,n,2 UW t when the degeneration is of pinching type and C 3n = C giiTlu i UW ( U C g2:Tl2i i with g x + 
gi = g and n\ + 712 = n, with 1 — 2^ — rii < 0, when the degeneration is dividing. The components 
Cg,n,h are nere Riemann surfaces with genus g, n punctures and h non overlapping disks removed 
which will become the punctures in the degeneration limit. The fact that the holomorphic integrals 
react correctly under the degeneration of the curve C 9 , n is a remnant of the analogous property of 
the conformal field theory ll28l and is indeed a consequence of the construction we performed in 
the previous section. We assume the shrinking cycle do not intersect the contour system along 
which (12.191 ) is evaluated. 
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Let us focus on the dividing case first. In this case the prime form E(z, z') behaves as follows. 
If both its arguments belong to a given same component, the prime form reduces to the prime 
form on that component, while if its arguments belong to different components, then E(z', z") ~ 
Ei(z\ P 1 )E 2 (P 2 , z")t~ l l 2 , where P\ 2 correspond to the punctures created by the dividing. To see 
what the prime form degeneration implies for the generalized matrix model measure and potential, 
we have to split the integration contours in components according to the dividing decomposition. 
This splits the {z;} in two sets according to which components of the contour they are integrated 
along, namely N' of them on the first component and N" on the second with N = N' + N" . 
Correspondingly, also the puncture set will split in two subsets w = w'Uw" , one for each component. 
By using the above degeneration formulas for the prime form and the fact that the holomorphic 
harmonic differentials cu a reduce to the ones relative to the two splitting factors, we get that 

Z^' n (w, m, p, v) ~ Z C N 3 r i+1 (w' U Pi, m' U ml, p', v')Z%?; n2+1 (w" U P 2 , m" U m* 2 , p", v") (3.1) 

where m* = bg s N'-J2 k , m k > +g s Q(gi - 1) and m 2 = bg s N" -J2k" m k"+9sQ(g2-l) after using 
momentum conservation. In the computation of (13.11) one needs to count the two extra zeros for the 
reference holomorphic differential to be placed at the location of the two resulting punctures. The 
direct computation of the above mass formulas from the integral (12.191) indeed gives, for example, 
m* = —bg s N" — g 2 g s Q + Ylk" m k" where the first term comes from the generalized Vandermonde, 
the second from the measure term in the square bracket and the third from the punctures. The 
computation of m 2 is identical. Notice that m\ + m 2 = —g s Q corresponding to the fact that the 
two Liouville insertions generated at the punctures are conjugated and therefore the masses of the 
two flavors at the two punctures p are equal. The conformal modulus which gets traded to the 
mass is the total filling fraction between the two integrals. Formula (13.11) is valid for each dividing 
degeneration such that g = g\ + g 2 , n = ri\ + n 2 with 1 — 2gi — rii < 0. 

In the pinching case, the prime form restricts to the one of the degenerate surface together 
with the holomorphic differentials with zero a 3 -cycle. The left over one scales as 2iriu g (z) ~ 
d z log §7j^ U P to 0(t) terms. In this limit, since to in the conformal factor has been chosen to be 
regular, two of its zeroes will be at the two punctures generated at the pinching node and by direct 
computation one gets 

(w, m, p, v) ~ Z C N 9 ~^ +2 (w U {P, P 2 }, m U {m% m*_}, p, v) (3.2) 

where p a are the momenta in the g — 1 left over handles and 

m *± = ~\- ± iPg> ( 3 - 3 ) 

with p g the momentum in the squeezed one. The two contributions to the above mass formulas 
arise respectively from the term in the square brackets in the measure and the second term in the 
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potential. Once again, the two masses at the generated punctures are Weyl conjugated m* + + m*_ = 
—g s Q. In the pinching case the conformal modulus which is traded for the mass is the momentum 
flowing in the squeezed handle. The formulas above are general and valid at finite N. 

In the following subsection we will discuss in detail the punctured torus case as an illustration. 
One could also study for example the punctured genus two case. This case is special since all 
degenerations reduce to punctured tori *. We discuss some aspects of the Seiberg-Witten geometry 
on genus two curves in the Appendix. 

3.1 Degenerations of punctured tori 

In this subsection, we concentrate on the pinching degeneration of a torus which leads to a sphere 
with two more punctures. Associated with the torus with n punctures we can consider a class of 
quiver gauge theories [|2l [3J whose particular weak coupling descriptions include a gauge theory 
with circular quiver. Specifying a particular weak coupling description corresponds to choosing a 
particular marking (a pants decomposition) of the Riemann surface. This gauge theory will reduce 
in the pinching of the torus to a linear quiver theory with n — 1 SU{2) gauge groups associated 
with a sphere with n + 2 punctures. In what follows, we verify that the generalized matrix model 
correctly reduces to the Penner type matrix model on the sphere ffT3l . 

Since on the torus the canonical bundle is trivial, the choice of a base point is not needed. The 
prime form is E(z, w) = ^jrji^y- and therefore the generalized matrix model (12.191) reduces to 

r N 

j i=i i<j 

up to ^-independent factors. The potential is 

n 

W{z) = J2 2m k^ge 1 {z-w k )+47rpz, (3.5) 

k=l 

and 

oo 

01 (z) = 2sin{7rz) Y[{1 - e 2niz q m ){l - e- 2mz q m ){l - q m ), (3.6) 

m=l 

* The dividing degeneration C 2 ,n — > Ci,n'+l UCin-n'+i generates two punctured tori. Indeed, the genus two prime 
form in such a degeneration reduces to the relevant ^-functions on the two tori since the period matrix at genus two 
becomes diagonal in the degeneration limit. In the pinching case C2, n — > Ci.n+2, the genus two ^-function entering 
the explicit expression of the prime from as E(z,w) = 7 ^ T ) contracts to the torus ^-function times a 

contribution from the off-diagonal term of the period matrix which cancels in the degeneration limit the square-roots 
of the abelian differentials appearing in the denominator of the prime form. 
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with q = e 2mT . The momentum conservation is given by 

n 

- J2^k + bg s N = 0. (3.7) 



k=l 



Also, the identification of the moduli of the torus and the gauge coupling constants of SU(2) gauge 
groups q k is as follows Il5ll20l: 



e 2ni( wi -w 2 ) = ft) e M» 2 -3) = ?2) ^ e 2 *- 1 ^ = g n _i, e 2niT = Y[q k , (3.8) 

k=l 

which, by fixing w n = 0, leads to 

e — q n -\ = r n-i) e — q n -2Qn-i = £n-2, ...,e — (ft . . . g n _x = ii- (J.yj 

The mass parameters and one of the Coulomb moduli correspond to vn k and p respectively. 

Let us consider the pinching degeneration of the torus. We take Imr — > oo which corresponds 
in the gauge theory to the decoupling limit of the n-th gauge group q n — > 0. To consider the 
behavior of the generalized matrix model in this limit, we first observe that the prime form reduces 
as 

= (dcr 1/2 (dcr i/2 (c-o, (3.io) 

where in the last line we have changed coordinates to £ = e 2mz and ( = e 2mw . It is straightforward 
to see that the Vandermonde determinant of (13.41) reduces to that of the usual /3-deformed matrix 
model. The potential also reduces to 



W(0 = }^ 2m * lQ g(£ " **) + 2 (~9sQ/2 - ip) log f , (3.11) 

fc=i 

where we have used ( 13.91) with t n = 1. Note that the first term corresponding to the momentum 
at £ = comes from the measure factor u(z) 1+b2 of the generalized matrix model (12. 19b . By the 
pinching, the punctures at £ = and oo are created. However, the latter disappeared from the 
potential, which thus reduces exactly to the Penner type matrix model lll"3~l[T9l t 

n-1 

W{z) = 2m k - h) + 2m n \og(z - 1), (3.12) 

k=0 



^The convention here is slightly different from the one in lfT3l [181 [191 . The momenta are related as 2m = ?Bdv- 
Our convention leads to the momentum conservation (13 , 1 3b . 
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with the momentum conservation 



n 



^2m k - moo + bg 8 N = g s Q, 

k=0 



(3.13) 



where is the momentum inserted at infinity. The relation between the parameters tk and the 
gauge couplings [19] (See also 0) is the same as the one defined in (13.91) with t = 0. It follows 
from (13.111) that m = —g s Q/2 — ip. 

Let us then analyze the momentum conservation under this degeneration. On one hand, in the 
original generalized matrix model, the conservation is described by (13.71) . On the other hand, in 
the Penner type one, the conservation is (13.131) . The momentum at infinity is then 



These values of the momenta m and are the ones which were already derived in the generic 
analysis of the previous section ( I3.3I ). Note that there is a slight difference between the momenta 
m and m^, which however disappears in the large N limit. 

The original generalized matrix model has n — 1 independent filling fractions v a . Recall that 
the overall J2 a v a — bg s N is constrained by the momentum conservation. Thus, by adding p, we 
have n independent parameters which are identified with the vevs of the vector multiplet scalars. 
The degeneration limit and the above argument mean that p in the potential is the vev of the n-th 
SU(2) vector multiplet scalar and some combinations of the filling fractions are the vevs of the 
other SU(2) scalars. We will give the precise identification in the large N limit in the next section. 

4 Spectral curve of the generalized matrix model 

In this section, we derive the spectral curve of the generalized matrix model (12.191) in the large N 
limit and show that it coincides with the Seiberg-Witten curve of the corresponding gauge theory. 

In the large N limit, the evaluation of (12.191) reduces to the calculation of the critical points. 
The condition for criticality is given by 



where the potential W(z) is defined in (12.171) and we have used the momentum conservation (12.61) . 
Then, the prepotential in the large N limit, defined as exp {jF/gf) = Z, is given by 




(4.1) 




(4.2) 
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where each "eigenvalue" Zi satisfies (14 .lb . 

It is natural to assume that the eigenvalues are distributed in line segments around the critical 
points of W(z), similarly to the usual matrix model. Indeed the second term in (14.21) reduces 
locally to the standard Coulomb gas potential. We denote the line segments as C a where a = 
1, • ■ • , n + 2g — 2. We assume that C a do not include the base point z* and the punctures Wk, 
at which the potential W(z) diverges. We denote by N a the number of eigenvalues on the line 
segment C a , where N a satisfies YT^=i~ 2 N a = N. 

Let us introduce the eigenvalue density current p(z) supported on {C a } and normalized as 
§ c p(z) = bg s N a = v a . Using the variables introduced above, the prepotential and the condition 
for criticality are written as 

T=\ p{z)W{z)-f [ p(z)p(z') log ^1 (4.3) 

dW(z) - 2 [ p(z')d z log ( = 0, (4.4) 

respectively. Here, z in (14.41) is on either of the line segment C a and the integral is defined as the 
principal integral. 

In order to solve the above condition (14.41) , we define the following one form, which is the 
generalization of the resolvent of the usual matrix model 

fl(2 ),^ c/M4l0g (|m). (4,) 

This "resolvent" is defined at generic points z on the Riemann surface contrary to the second term 
in (14.41) . Note that the resolvent as well as dW(z) are single-valued one-forms on the Riemann 
surface. The resolvent has cuts at the line segments C a and a simple pole at z*. Also, the filling 
fractions are obtained by integrating the resolvent along the cuts as 

u a = <f R(z). (4.6) 
2^ J Ca 

On the line segments C a , the resolvent behaves as 

R(z + iee iv{z) ) + Rtz - iee^ {z) ) = 2P f p(z')d z log f ffi' Z '\ ) = dW(z), (4.7) 

JZ a c a \E{z,z*)J 

R{z + iee iLp{z) ) - R(z - iee lip{z) ) = ^p(z')d z log f |^2 j = -2irip(z) (4.8) 

where we take the real number e infinitely small and we assume that a properly defined function 
(p(z) exists such that z + iee lif ^ or z — iee ltp ( z > does not go across the cuts C a when z moves along 
C a - The integral in (14.71) is principal integration, which is given as an average of integral along the 
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path above the singularity and that below the singularity. The resolvent should be determined such 
that (14.71) and (14.81) are satisfied for z G C a . A candidate of the solution for (14.71) is 

R (z) = l -dW{z). (4.9) 

However, it does not reproduce the correct structure of singularity expressed in (14.81) . We need 
singular contributions: 

R(z) = ^dW(z) + R(z) sing , (4.10) 

where (14.71 ) and (14.81 ) impose 

R(z + iee^) sing + R(z-iee^) siQg = 0. (4.11) 
R(z + iee ilp{z) ) sing - R(z - iee ilp{z) ) sing = -2mp(z). (4.12) 

The above discussion is valid for a generic potential W(z). In the following, we use its explicit 
form (12.171) to determine the resolvent R(z). Then, we find that 

/ n g 
p(z')d z \ogE(z,z') - y^m fc rf z log E(z,w k ) - 2tt y^ j p a uj a (z) (4.13) 

C " k = l 0=1 

does not depend on the base point z*, where we used the momentum conservation (12.61) and ignored 
the subleading term in the large iV expansion. We see that R(z) sing has cuts in the regions C a and 
simple poles with residues nik at z = Wk- Moreover, it is independent on the base point z* as 
expected. 

From (14.111) . we see that the sign of -R(^) S in g changes across the cuts. Therefore its square only 
displays singularities at the punctures z = w^. From (14.131) we see that these are at most quadratic 
poles with coefficients m^ 2 . The spectral curve of the generalized matrix model thus reads 

n 

Rsm g (z) 2 = J2 m kMz,w k ) + C(z) (4.14) 

k=l 

where r](z, Wk) are quadratic Strebel differentials, with double pole at uik, and ((z) is a quadratic 
differential which has at most simple poles at w k . ((z) is determined in terms of n + 3g — 3 
parameters; in particular it depends on the n + 2g — 3 independent filling fractions u a and the 
g internal momenta p a . This form of the curve is the same as that of the Seiberg-Witten curve 
of quiver gauge theory as a cover of the base Riemann surface BUO. As discussed in ||29l , the 
physical information is included in the Prym variety of the Seiberg-Witten curve, rather than the 
Jacobian variety. This reduces the number of independent periods of the Seiberg-Witten differential 
to n + 3g — 3, which agrees with the number of the filling fractions and p a . 
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In order for the above argument to be a check of AGT correspondence at all genera, we need 
to show that the spectral curve is indeed the same as the one proposed, in the Virasoro conformal 
block side, to be the Seiberg-Witten curve: the insertion of the energy-momentum tensor in the 
conformal block 



T(z)l[V ak (w k )) (f[Va h {w h )) (4.15) 




in the semi-classical limit, which corresponds to the large N limit in the generalized matrix model. 
It was claimed that xdz coincides with the Seiberg-Witten differential. This was checked in the 
cases on a sphere and a torus already in [|5]|. A useful way to capture the energy-momentum tensor 

<jy(z) 

insertion is to consider the insertion of the degenerate field $1,2(2) = e fa in the conformal 
block. By this insertion, the conformal block satisfies the BPZ equation and can be expanded in g s 

as m 

^ 2 (4^w ~ exp (5 + ^r/ x ^ dz>+ --) (4 - 16) 

fc— 1 / Cg 

Thus, the counterpart of the Seiberg-Witten differential in the Virasoro side can be obtained by 
calculating 

a. (<M*)nLi^>*)) c , 

x{z) = bg s — log — - (4.17) 

dz (Uk=i v «k( w k))c a 

in the semi-classical limit. 

In the following, we rewrite (14.171) in terms of the generalized matrix model by using the 
discussion in section |2] The corresponding calculation in the case of the sphere has been done in 
||2TT| . We use the holomorphic half of the integrand in (12.201) to evaluate the conformal block in 
(14.171) . The conformal block with the insertion is also obtained just by changing n to n + 1 and by 
regarding that m n+ \ = % and w n+ \ = z. The momentum conservation (12.61) is slightly modified 
to 

S mfc + Yb + 9sQ{1 ~ g) = b9sN - (4 ' 18) 
fc=i 

By collecting the factor dependent on z, we obtain 

8 , / , , 1 -1— r _,/ N Q . , . at Sfc m fe _WM_ 

bg S Tj~ log [ u(z)^ 1 [[E(z, ^)^E(z, z) b as e ^ 

OiN- 1 ), (4.19) 
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where we have used the deformed momentum conservation (14.181 ) in the second equality. The 
expectation value is defined by 

(...) = — — s . (4.20) 

A subtlety is that the numerator in (14.201) is defined with the deformed momentum conservation 
(14.181) while the denominator with the original momentum conservation. However, note that the 
effect of the insertion of the degenerate field is a subleading contribution in the large N limit, as 
in (14.161) . Note also that the deformation of the momentum conservation gives rise to that of the 
prepotential Fq, which does not depend on z and disappears by the partial derivative in terms of 
z. Thus, we can use the same external momenta m k and the same iV as those of the denominator 
of (14.201) to evaluate the numerator in the expectation value. Also, some of the internal momenta 
in the numerator are shifted by ±6/2 depending on where we insert the degenerate operator as 
discussed in iU, but again, the effect of this shift does not produce the factor dependent on z in the 
large N limit. Thus, the expectation value in (14.191 ) can be evaluated by substituting the solution 
of the condition of criticality (14.11) or (14.41) . which leads to 

From the discussion above, we finally obtain 

x(z)dz ~ -^M + R ( z ) = J R sing (^), (4.22) 

up to 0(N~ l ) terms. This shows that the resolvent -R S in g is indeed x{z)dz in (14.151) . which was 
claimed to agree with the Seiberg-Witten differential. 



4.1 Dependence on internal momenta 

In this subsection, we show that the spectral curve of the generalized matrix model depends on 
the Coulomb moduli parameters in the same way as the Seiberg-Witten curve, which completes 
the check of AGT correspondence at all genera. In order to do that, we first have to specify the 
marking of the Riemann surface C g>n , which is done by choosing n + 3g — 3 physically independent 
cycles. This marking determines the conformal block labeled by a trivalent graph on one hand, and 
also the corresponding weak coupling description of the gauge theory on the other hand. On the 
gauge theory side, each vacuum expectation value of each vector multiplet scalar is obtained by an 
integral along each cycle of the marking 7 r : 

a r = w- { Asw, (4.23) 
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where A$w is the Seiberg-Witten differential. We would like to show that the spectral curve also 
satisfies the corresponding relation. 

Before going to that issue, we show in the following discussion that each pair of pants has one 
cut under the marking. In general, three-punctured sphere always has one cut lfT5ll because the 
classical potential of the corresponding Penner type matrix model with generic mass parameters 
has one extremum. We can also show this by considering the dividing limit of four-punctured 
sphere: under the decomposition, it is known that one cut exists in the pants including the punctures 
originally at and q, while the other cut exist in the pants including the punctures at 1 and oo. 
Furthermore, unless the mass parameters vanish, there are no punctures belonging to the cut. This 
indicates that three-punctured sphere with generic mass parameters has one cut apart from the 
three punctures. Since we have shown that our generalized matrix model behaves correctly in the 
degeneration limit, the cuts of its resolvent should be placed in such a way that each pair of pants 
has one cut. 

Once the marking is specified, an explicit correspondence should be determined between the 
n + 3g — 3 internal momenta of the Virasoro conformal block and the n + 2g — 3 independent 
filling fractions together with the n parameters p a in the potential. In order to see this explicitly, 
we can use the results on the degeneration obtained in section[3] Note that among the intermediate 
states in the conformal block, only the primary state remains and appears as the external state in the 
degeneration limit. Thus, the external momentum appeared by the degeneration which shrinks a 
particular cycle can be seen as the internal momentum of the corresponding place before the degen- 
eration. In order to determine all the internal momenta, we pinch all the A-cycles of the Riemann 
surface C 3i „ until it becomes Co, n +2# and then divide Co, n +2g until they split into three-punctured 
spheres in a way that it reproduces the specified pants decomposition. Whichever marking we take, 
we can read off all the internal momenta by considering the corresponding degeneration limit. 

The A-cycles of C g , n can be chosen to coincide with some of 7/s. Thus, the set {7,.} can be 
divided as {7,.} = {A a } U {7 a } where a = 1, . . . , g and a = 1, . . . , n + 2g — 3, such that bi- 
cycles of C 9t n are A a and 7 a are the remaining cycles. By shrinking v4 a -cycles, we obtain (13.31) . 
Especially, in the large N limit, the internal momenta (multiplied by g s ) corresponding to the A a 
cycles of the Riemann surface C g are given by 

a a = ip a - (a = 1, • • • , g) (4.24) 

After pinching all the A a cycles, the remaining internal momenta can be obtained as written just 
below (13.11) by further considering the dividing limit which shrinks 7 Q cycles. Since each pair of 
pants has one cut, on which Np eigenvalues {zi} exist, the number of eigenvalues N' or N" on 
each side of the dividing component of the Riemann surface can be written in terms of sum of the 
filling fractions. In the large N limit, the internal momenta (multiplied by g s ) flowing through 7 Q 
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cycles are given by 

a « = - ^2 m k + ^(±ip a ), {a = l,--- ,n + 2g-3) (4.25) 

/3 fc a 

where the range of the sum depends which "f a cycle we shrink under the marking of the Riemann 
surface. The third term is the contribution from new punctures which appear by pinching A a cycles 
in the previous step. The sign ± is determined for each puncture according to the direction of the 
A a cycle. Although the internal momenta (14.241) and ( 14.251 ) are obtained in the degeneration limit, 
we assume that these relations hold for arbitrary moduli parameters of the Riemann surface C g , n 
because the filling fraction v a and the parameter p a are independent of the moduli parameters. It 
is remarkable that the form of our generalized matrix model in the large N limit is universal for 
any choice of markings of the Riemann surface. However, the identification of the parameters with 
internal momenta depend on such choice, which reflects in the difference of the conformal blocks 
labeled by different trivalent graphs. 

So far, we have discussed the relation between the internal momenta of the Virasoro conformal 
block and the independent filling fractions and parameters p a . In the following, we check that the 
corresponding relation as (14.231) is reproduced from the spectral curve under these identification of 
the parameters. By using (14.131) , we can explicitly calculate the A a -cycle integrals as 

-Rsing = ip a (a = 1, • • • , g) (4.26) 

Note that the convention here is § A uj}, — 5%. Taking into account the identification of the parame- 
ters (14.241 ), we see that the corresponding relation as (14.231 ) for A a -cycle is confirmed. 

For the 7 a cycles, the check is not direct, because 7 a are not the contours C a along the cuts of 
the spectral curve of the matrix model. Recall that there are indeed n + 2g — 2 cuts in the spectral 
curve, but the integral of the resolvent along one of them is not independent due to the momentum 
conservation. We will call the remaining n + 2g — 3 cuts the independent ones. It follows that the 
cycles around the independent cuts, which we have denoted by C a (a = 1, • • • n + 2g — 3), can 
be expressed as sum of the three cycles which go around each leg of the pants. That is, C a can be 
written as a linear combination of 7 r (r = 1, ■ ■ ■ n + 3g — 3) and D k (k = 1, • • • n), where D k are 
the contours around the simple poles of R sing . Thus, 7« (a = 1, • • • n + 2g — 3) can be written as 
7a = n a a A a + c^Cfi + s a h Dk, where the coefficients n a a , c a @ and s a k are integers. (Note that 
in addition to them there could exist a cycle around the base point. However, when we consider 
_R sing , this dependence disappears in the large N limit, as stated in (I4.13I ).) In terms of these, we 
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7n-l C fl 

Figure 1: A marking of Ci n corresponding to the elliptic quiver and the A\- and 7 Q -cycles. E.g., 71 = 
Ai + Ci + £>i and 72 = Ai + Ci + C 2 + £>i + D 2 . 



calculate the integral over 7^ cycles as 

R(z) -^m k u Wk)Z *{z) - 2n^2p a uj a ( 



— I R si ng(z) = — 

2m L 2ni , 

^ a J la l k= l a = l 

= ^c a p iJp-^s a k m k + i^n a a p a , (4.27) 

k a 

where u X:y (z) = d z log ^f^y is the Abelian differential of the third kind. 

To determine the coefficients, we briefly show the example below. In the case of the torus C\ i7l , 
there is only one A-cycle and the number of independent cuts is n — 1. For n — 1, there is no 
independent cut and the parameter p is identified with the vacuum expectation value of the vector 
multiplet scalar of J\f = 2* gauge theory as analyzed in [|20l . For n > 1, we consider the marking 
corresponding to SU (2) n elliptic quiver gauge theory and choose the Ai-cycle and 7 a -cycles as 
depicted in Fig \T\ As stated above, each pair of pants has one cut. (This can also be explicitly 
checked by solving the equations of motion in the weak coupling limit.) Therefore, it is easy to 
obtain the coefficients above nl = 1 and 

s k a = 

I U [p = a + 1, . . . ,n - 1J 

This indicates 





1) a (Jfe = a+l,...,ra-l) 



1 f a a 



T« /3=1 k=l 



in addition to the one obtained by the A-cycle integral (14.261) . Comparing it with the result obtained 
from (14.251) . we see that it exactly reproduces the internal momenta flowing through 7^ cycle. 
Similarly, it is straightforward to check for generic case that (14.271) coincides with that obtained 
from (14.251) . Thus, we have confirmed that the corresponding relations as (14.231) for 7^ cycles are 
reproduced. 
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5 Conclusions 



In this paper, we have shown that the perturbative analysis of Liouville correlation functions dis- 
plays in the large N limit holomorphic factorization of the surface integrals and leads to generalized 
matrix models, defined on the cover of C g>n , which describe the relevant Virasoro conformal blocks. 
We provided an all genera check of the AGT correspondence by obtaining the Seiberg-Witten data 
from the saddle point analysis of these generalized matrix models. 

We underline that the models presented in this paper could be useful for the exploration of the 
full set of gauge theories with generalized quiver structure of 01. Indeed, so far most of the analysis 
of the AGT correspondence has been focused on the linear and elliptic quiver cases, mainly due to 
the lack of calculational tools for higher genera. However, to fully exploit the generalized matrix 
model approach one should be able to extend its analysis to finite N. This would amount to provide 
a full derivation at finite iV of the holomorphic factorization, which in turn would give a precise 
prescription for the contour integrals possibly extending the recipe of |[T6ll22l to higher genera. 

Notice anyway that our large N analysis depends only on the homotopy class of the contours 
via the filling fractions as explained in Section 4.1. In this sense our results are universal with 
respect to a particular choice of contour's representatives in the matrix integral. 

Moreover, for gauge groups of higher rank, which according to lfT3l [T4l should correspond to 
multi-matrix models, this approach could shed light on the description of strongly coupled sectors 
naturally appearing in the general framework and not admitting a known lagrangian description. 

Another very interesting issue to explore is the relation of the generalized matrix models with 
the quantization of integrable systems ll30l[TTT| . (See also OTI .) In particular this could provide an 
alternative derivation of the quantum Hamiltonians for Hitchin integrable systems and generalize 
it to higher genera. 
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Appendix 



A Addendum: playing with genus 2 curves 

A hyperelliptic curve C g of genus g is given by the equation 

V 2 = P2g+2(x) 

where P2 g +2 is a polynomial of degree 2g + 2 and is realized * in the total space of TP^ . 

As it is well known, all genus 2 curves are hyperelliptic. These are realized in general by a 
sextic polynomial equation 

6 

i=i 

which we denote by C 2 - 

The complex structure moduli M.2,0 of genus two curves is then obtained by considering the 
complex parameters {a,;} modulo the action of the permutation group S 6 and the PSX(2, C). 

A basis of abelian differentials is given by u a = xa ^ dx , a = 1,2, while a basis of quadratic 
differentials is given by Q = a = 1.2,3 11321 

The Seiberg-Witten (SW) geometry of the SU(2) theory at genus 2 is specified by a double 
cover of C2 in T*C2- As such, this is specified by a general quadratic differential on C2 in the form 

w 2 = $ 2 (A.2) 

where $2 = J2 a Ka4>a can be expanded in the Coulomb moduli K a of the theory. 

The perturbative expansions of the theory are available in the vicinity of the degeneration locus 
of the moduli space, namely around 

^2,0 = ^1,1X^1,^^,2 (A.3) 

The second factor in (IA.3I ) is still generically not lagrangian and has to be degenerated as dA4i j2 = 
.M^i x jMo,3 U .Mo,4 to reach corners around which known lagrangian descriptions are available. 
The first factor in (IA.3I) . being given by two copies of the J\f = 2* 577(2) theory, is already 
lagrangian. The first degeneration is dividing and the second one is pinching. 

Let us discuss the dividing case in detail. This is reached by taking the limit in which three 
branch points in (IA.1I) collide. To be concrete, let's fix the position of two of them at and 00, 

■^Under conformal inversion x = 1/x' on the Riemann sphere, the stability of the description is guaranteed by the 
transformation y = y'x' _ ' 9+1 ' = y' (§§7) 2 ■ 
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write our curve as 

y 2 = x(x — a\e)(x — a 2 e)(x — a 3 )(x — a 4 ) (A.4) 
and take the limit as e — > 0. The curve (IA.4I) becomes in the x coordinate 

y 2 = x 3 (x - a 3 )(x - a 4 ) 

which, redefining y = xy, reads 

y 2 = x(x — a 3 ){x — a 4 ) (A.5) 

that is the torus with a puncture at x = 0. 

Let us check now the SW geometry in the degeneration limit. The issue to discuss is just 
the scaling of the Coulomb parameters in this limit. Notice that the degeneration is obtained by 
contracting to zero two nearby branching points, therefore we are saturating a complex structure 
modulus corresponding to a Beltrami differential fi e with support around the origin of size ~ e. 
This is dual to a holomorphic quadratic differential which, not to have a vanishing overlap integral 
with ji e should not be zero at x = 0. This is uniquely determined to be Therefore, along the 
limit with e — > 0, the corresponding parameter in the SW curve has to scale away. 

As a consequence of the above reasoning, exposing the e-parameter, the SW curve is parametrized 

as 

2 u'e + m 2 x + ux 2 . , , 2 / A ^ x 

w 2 = (dx) 2 (A.6) 

y 2 

The degeneration of C 2 is easily kept into account in the SW geometry which becomes 

,2^-1 

J 2 

The standard parameterization of the punctured torus is in the coordinates where the puncture sits 

at oo. Therefore, we rewrite the elliptic curve (IA.5I) after the inversion x = to pull the puncture 

i 

W? 



9 m x + u , , n9 
w 2 = zj^(cfe) 2 (A.7) 



at x' = oo and redefine accordingly y = -A^y'. After this, the SW curve reads 



which we can put in the representation with respect to the periodic coordinate via the Weierstrass 
parameterization 5 x' = V(z) + c and y' = ^V(z) so that we stay with 

w 2 = [m 2 V{z) + (u + cm 2 )] (dz) 2 (A.8) 

which is the SW curve for a copy of the J\f = 2* theory. 



s The constant c needs to bring iA.5i to the standard Weierstrass form where the quadratic term vanishes and can 
be computed explicitly. 
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The other copy corresponds to the other half in which the original genus 2 surface was split. 
Let's see how to get this second copy. In order to do it we have to consider the curve in the coordi- 
nate appropriate for the other half, namely we have to change (IA.4I) to x = e/x and correspondingly 

3 /2 

V = Vhw a f ter which we get 

\ x ) 

y 2 = x(l — a\x)(l — a 2 x)(e — a 3 x)(e — a±x) . (A.9) 
(IA.9I) becomes after the degeneration limit the curve 

y 2 = x 3 (l — a\x)(l — a2x){a^ai) 
which we bring to the form of a punctured torus by redefining y = yx and get 

~2 

y = x(l — aix)(l — a2x){aza$) . 
Let's follow now what happens to the SW curve (IA.6I) in the x patch. This becomes 



2 u' + m 2 x l +uex 2 2 
w = ^ (dx) 2 (A.10) 

y 2 

which in the limit e — > has the same form of (IA.7I) . but referring to the second punctured torus 
with an independent Coulomb parameter u'. So, following the same procedure leading to (IA.8I) . 
we get the second copy of M = 2* with an independent gauge coupling and Coulomb parameter 
but the same mass as the first. 

A pinching of the genus 2 curves (|A.4b can be obtained by letting e — > — for example. In such 
a case the curve gets to 

y 2 = x(x — aia3/a 2 )(x — a 3 ) 2 (x — a 4 ) (A. 11) 

which is, after renaming y = y(x — a 3 ), the twice punctured torus 

y 2 = x(x — aia 3 /a2)(x — a^). (A. 12) 

Correspondingly, the holomorphic quadratic differential entering the Seiberg-Witten curve 
(IA.6I) becomes 

u'ari/ao + m 2 x + ux 2 , , n9 
$2 ^2 ( dx ) (A. 13) 

which explicitly displays quadratic poles at the two images of x = a 3 with equal coefficients. 
This coefficient is actually fixed by the Coulomb parameter corresponding to the ungauged group 
SU(2) at the end of the shrinking of the handle. 
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